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£h ! ABSTRACT 

We consider the quantum Friedmann equations which include one-loop vacuum fluctuations 
due to gravitons and scalar field matter in a FLRW background with constant e = —H/H 2 . 
After several field redefinitions, to remove the mixing between the gravitational and matter 
^+ ! degrees of freedom, we can construct the one loop correction to the Friedmann equations. 

Due to cosmological particle creation, the propagators needed in such a calculation are 
typically infrared divergent. In this paper we construct the graviton and matter propagators, 
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making use of the recent construction of the infrared finite scalar propagators calculated on 
a compact spatial manifold in [lj. The resulting correction to the Friedman equations is 
suppressed with respect to the tree level contribution by a factor of H 2 /mi and shows no 
secular growth. 
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1 Introduction 

Since the classic work of 't Hooft and Veltman on quantum gravity on flat space-time [2], 
many authors have studied the quantum behavior of gravitons. In particular, because of 
the potential relevance for inflationary cosmology, the quantum behavior of gravitons on a 
(locally) de Sitter background has been a widely studied subject over the past years [3J HJ EJ 
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de Sitter space is interesting for its relevance for inflationary cosmology. Therefore there are 
several works that study the potential influence of quantum behavior on inflationary observ- 
ables [301 [3H [32j [33] . Another line of research deals with the back-reaction of gravitational 
waves on the background space-time [THl 03 [T21 [20] ■ However of more interest for the present 
work is the one loop back-reaction by virtual gravitons on a de Sitter background, which 
has been calculated by several authors, using different techniques [T3l [28J, [29]. Since it is not 
clear whether in these works exactly the same quantity is calculated and the renormalization 
schemes differ, the numerical coefficients differ. However the main result is the same: one 
loop graviton contributions to the expectation value of the energy momentum tensor result 
in a finite, time independent shift of the effective cosmological constant. Since the contri- 
bution can always be absorbed in a counterterm [28] , the exact numerical coefficient coming 
from such a calculation has no physical meaning. 

The goal of this paper is to go beyond the works mentioned above and calculate the one 
loop effective action induced by gravitons in a more general background space-time using 
dimensional regularization. The geometry we consider is a Friedmann-Lemaitre-Roberston- 
Walker (FLRW) geometry with Hubble parameter H = - and the additional constraint 
that 

—I « 

is a constant. If the universe contains only one type of matter, e.g. dust, radiation or a 
cosmological constant, this constraint is satisfied. In particular in matter era we have that 
e = 3/2, in radiation era e = 2 and de Sitter space corresponds to the limit e — > 0. One 
immediate problem with working in such a space-time, instead of de Sitter, is that for con- 
sistency of the Einstein equations the addition of matter fields is unavoidable. Whereas in 
de Sitter space, the only relevant metric fluctuations are the tensor modes (gravitational 
waves), in a more general setting one has to take the mixing of gravitational and matter 
degrees of freedom into account [311 EH [231 I2U [251 [26] . 

Naively one might think that any back-reaction will be insignificant, since it will in a 
cosmological setting typically be suppressed by H 2 /m 2 with m p the Planck mass. The reason 



that this is not necessarily the case has to do with a second complication. This complication 
is that gravitons on a cosmological background space-time behave, apart from the tensorial 
structure, somewhat similar to massless scalars, with a certain amount of coupling to the 
Ricci scalar. Such a scalar field however is long known to posses problems in the infrared [361 
[37] . for a wide class of cosmological backgrounds. What happens is that due to cosmological 
particle production correlations of the Bunch-Davies vacuum grow too fast, causing the 
propagator to diverge in the infrared. Such a divergence should of course be regulated and 
in a recent work [1] the massless scalar propagator on any constant e cosmological space-time 
was constructed, by assuming that the universe is described by a spatially compact manifold. 
Effectively this simply implies that any momentum integral has an infrared cut-off at some 
scale &o- in this work, we shall use the propagators constructed in [TJ to obtain an infrared 
finite answer. 

While the answers obtained in this way are infrared finite, particle production is of course 
a physical phenomenon. Long range correlations will therefore still be enhanced as time 
goes on. Therefore initially small quantum fluctuations, might - at least in principle - grow 
significantly in time. The contribution of these fluctuations to the stress energy tensor might 
therefore also grow in time, leading to a potentially significant back-reaction via the Einstein 
equations on the background space-time. 

It is precisely this effect, that has led to several papers concerning quantum contribu- 
tions to the stress energy tensor in the case of de Sitter space-time [61 [101 [H [23j [2U [251 138] . 
In these works it is found that at two loop order, the backreaction due to gravitons might 
become significantly, while for scalars this might happen at three loop order. Such behavior, 
if indeed physically viable, has profound implications for the cosmological constant problem 
(see e.g. Ref. [39] for a review). In this paper we will not consider these higher loop effects, 
but we shall perform a one loop calculation, but in a more general space-time. 

In section [2] we briefly review our background geometry. In section [3] we briefly summarize 
the construction of the propagator as given in [TJ. In section H] we apply these results to 
calculate the relevant kinetic operators and propagators in a model with both matter and 
gravitational degrees of freedom. In section 14.21 we calculate the one-loop effective action 
contribution to the quantum corrected Friedmann equations and renormalize the theory. We 
discuss our results in section 14.31 and we conclude in section 



2 Geometry 

The background space-time we consider is is the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) geometry in conformal coordinates 

9tw = a 2 V^ ; Via> = diag(-l, 1, 1, 1), (2) 

with the additional constraint, 

H TT a . . 

e = = constant ; rl = -. (3) 

H 2 a 

Here a dot indicates a derivative with respect to cosmological time t, related to the conformal 
time rj by dt = adrj. The FLRW geometry obeys the Friedmann equations 

3H 2 1 2H 1. 

~Pm = ; -(Pm+Pm) = 0, K = l6nG N , (4) 

K I K I 

with Gn being the Newton constant, pm and pm are the energy density and pressure due to 
matter. If one writes 

Pm = wmPm, (5) 

one immediately finds that ([3]) implies that u>m is constant. One can solve (jlj) for a to find 



a(V) = \(e - l)H 7]J ; e = -(l + w M 

(, , W(l-e) 1 ' 

H = H [{e-l)H oV ) =H a- e ; a V = -j^J 



(6) 



Notice that if e < 1, rj is negative and the expansion of the universe is accelerating. If e > 1, 
7] is positive and the expansion is decelerating. Hq is chosen such that the e — > limit of 
H corresponds to the one given in [22]. Given two points x and x, the following distance 
function will prove to be useful 

y(x; xj = = — {—(\V — m—ie) + \\x — x\\ ) . (7) 

7777 ?7?7 

Here the infinitesimal e > refers to the Feynman (time-ordered) pole prescription. In the 
following we shall omit the explicit argument of y(x;x), since this will be clear from the 
context. In de Sitter space y is related to the geodesic distance I as y = 4sin 2 (|iJ/). If 
y < 0, points x are timelike related to x, and if y > 0, they are spacelike related. We 
define the antipodal point x of x by the map rj —>■ —rj [10] ■ Notice that, since as long as e is 
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Figure 1: The causal structure in the conformal coordinates (J2J). The plot assumes e < 1, so 
the coordinates cover only the region rj < 0. The wavy line at rj = indicates future infinity. 
The lightcone of the point x is given by y = 0. If y — 4, the point x lies on the light cone of 
an unobservable image charge at the antipodal point x. 



constant, rj is either positive or negative, this point is not covered by our coordinate patch. 
If y = 4, x lies on the lightcone of an (unobservable) image charge at the antipodal point x, 
see figure [TJ 

The Levi-Civita connection is, 



r^ = a - [sy? + 5 a X + 8%vr] , 

(J/ 

while the curvature tensors are given by 

ir „ ,„ ( - - 2 ( - ] ) fej<sg - 5%5% p - sffisl + qssv,* 



fj,/3u 



a \aJ 

I- 2 



\\w 



R„„ =( ^ - 2(^) 2 ) („,„ - (D - 2)^°) + (^) 2 (Z) - 1)„, 



•-a^n^-D+r^v-D 



i2 

.o° Va z / / Va z 

where a' = da/drj and D denotes the number of space-time dimensions. 
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3 Scalar Propagator 

Since it will turn out that the graviton propagator can be expressed in terms of massless 
scalar propagators with different amount of conformal coupling, we need to find a solution 
for the following scalar Klein-Gordon equation 



-g (□ — £R) tA(x; x) = iS (x — x 



(10) 



where £ is a dimensionless scalar field coupling to the Ricci curvature scalar, and □ denotes 
the d'Alembertian, which reads when acting on a scalar 



If e is constant we can rewrite this equation as follows 



rTd^ + ^lv 



rr 



1 



(aa) 2 l iA(x;x) =i5 D {x — x) 



(11) 



(12) 



where a = a(x) and 



v 2 =(^ Y ±y-(l-e)- 2 ((D-l)(D-2e)Z- 1 -(D-l)(D-2)e+j(D-2)e 2 ). (13) 

The properties of such a scalar field propagator on a de Sitter background (e = in our 
notation) have been studied in the past [301 EJ 132] . 

This propagator in a general constant e space is given in pQ , by integrating over the mode 
functions and we simply quote the result here. The massless scalar propagator is infrared 
divergent if the index v (which is defined in ({TBI) is less then [D — l)/2. Therefore we work 
on a compact spatial manifold, with comoving radius k^ 1 . In this case, the integral over the 
mode functions becomes a discrete sum. However, if k is small enough, it is generally valid 
to approximate this sum by an integral, with an infrared cutoff. This approach has passed 
sever consistency checks on de Sitter space [61 H3J S3]- In the final propagator we can then 
recognize two different types of contributions: one coming from the 'infinite volume' integral 
(thus the integral, without an infrared cut-off), and one from the subtraction of the lower 
range of the integral. 

oo oo 

iA(x; x) = zA 00 (a;; x) + ^ 5iA N (x; x) + ^ MA N (x\ x) . (14) 

N=0 N=0 

Here the infinite volume propagator is 

, , "., , , ( D-\ D-\ D y 

**~tez) = ^ ^ 2 i^— +„,—-„; T ;l-- 

(15) 
and the correction terms due to the infrared divergence are 



\l-e) 2 HH 


r(^i- 


^)r(^i- 


-v) 



5iAn[x\ x) 



2—1 
HHCi - -C) 2 ) " 2z 2N+(D-l)-2u 



"^ fr ^2 k/ ^2i-N ( 16 ) 

X wi ' 'Jn-u > > Owl ~ 



(4tt)t 2iV + (£> - 1) - 2i/ 

8iA N (x;x) = 5iA N (x\x)\ u ^- v , 
where z Q = k Q \r)\ and 

i\N i r(^)p(i-z/ 

4J fc!^!(iv-jfc-^)! r(fc+^i)r(^+i-i/)r(iv-jfc-^ i-v) 



°M=(-7) ,,^,r , .MX,,,.. P-l^,, ..^ ,. — ~ • ( 17 ) 
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Notice that in an accelerating space-time (e < 1) r\ approaches zero at late times, while in 
a decelerating space-time (e > 1), 77 approaches infinity. This is an interesting observation, 
since it implies that in an accelerating universe, iA^ contains terms that grow in time 
if v < ^y~- In a decelerating universe both z'Ajv and iA N can contain growing terms. 
It is precisely these growing terms that can lead to a secular growth of quantum effects, 
described in the introduction. There are also other regularizations possible for the infrared. 
For example in [35] the infrared is regulated by matching an infrared finite space-time to 
the space-time under consideration. This ensures that no infrared divergences can occur. 
Results calculated using this regularization are similar in accelerating space-times, but show 
differences in a decelerating space-time. 

4 The gravitational and matter propagators 

For our model we shall consider the action of gravity plus a scalar field <f) = <p(x) with an 
arbitrary potential V(<p) 

^/^(^|^-iW) 2 -^)), (18) 

where k = IQtcGn = Wir/m 2 ? denotes the (rescaled) Newton constant, mp ~ 1.2 x 10 19 GeV 
is the Planck mass, A denotes the cosmo logical constant, D is the number of space-time 
dimensions and (<90) 2 = g^ u (d^(j))(d u (j)). By an appropriate choice of the potential V, such a 
model can mimic any mixture of fluids which are relevant for the evolution of the univese. 
We consider a conformally flat FLRW background and split the fields in a background 
contribution and a quantum contribution [7j[26j 



y flVJ 



gv» = g ^( v ) + Sg" u = a- 2 {rf v - ^a 4 ^) + 0{^ 2 ) (19) 

= $(??) + 0, 

where Sg^ = h^ u denotes the graviton field, ip^ u is the pseudo-graviton field and 5g^ u = 
—h^ u + h' J a h otu + 0(h 3 ). Notice that indices on the pseudo-graviton field ?/v are ra i se d and 
lowered with the full background metric g^uiv) = c^Vv- The background scalar field $ is 
homogeneous and thus only depends on (conformal) time. The background fields obey the 



tree level Friedmann equations and the scalar field equation of motion: 

H * ~ ^T A - ( D - in - 2) (2^ + y W) = ° 
a -i tf + £zi fl ._i A + _JL_(^ a _v W )=o (20) 

fly 

$" + (£> - 2)a#$' + a 2 — -($) = , 

from which one can derive the following identities 

y/i&' = y/2(D - 2)eaH ; \/^>" = ^2(D-2)e(l - e)a 2 H 2 + O(e') 

V^(*) = -V2(D-2)e(D - 1 - e)H 2 ; |^($) = 2(D - 1 - e)etf 2 + O(e'). 

(21) 

Notice that in slow roll inflation e' is non zero. In order to facilitate a comparison, we give 
here the slow roll expressions for the second and fourth relation (the first and third do not 
change) 






(22) 



dcp 2 " ' 2 

in terms of the slow roll parameters 



- 1 f 9V Y __2d 2 V 

ev= K\Vdj) ; Vv = ^Vdf 2 ' ( } 

Notice that in the slow roll approximation e and ey are identical. 

For the purpose of this paper we are only interested in the quadratic perturbations in 
the fields. After many partial integrations we find u 

C (2) =a D + ^ I ( Qs _ yy) (\ 6 P i6: _ l^rr) + *<?« - yri^rTJ V> 

- a D ~ 2 ^ (v^$") </> + tP-*Jhr^Z* (24) 

1 n,/„ a 2 ^^, _„„ \ , 1 



+ -a D <p(n s ~ w m + a-V K )<P+ -V^s^FaFp , 



1 An analogous result can be found in Ref. [26]. The main difference is that our result (f24l) includes also 
terms that vanish on-shell. 



where we defined 



F a = a 2 V M (< - itf&TVw) - <&V*% 



W = 2(D-2) 



a-'H' + ^-^H 2 - -A 



K 



_$' 2 _ v($) 
2 2{D-2)\2a 2 K } 



X 

y 



-2/R/2 



-na $ 



-lrj/^ 



-(H 2 -a- L H 



D — 2 



2 V a</> 



and 



□., 



"5 



-.d.y/^g^d, 



is the d'Alembertian as it acts on a scalar field. We add a gauge fixing term 



£gf = -^V-99"' r a r p 



and therefore we also need to add a ghost lagrangian 

Cw = -\f ZI ga~ 2 V lx 8F IJi , 



(25) 



(26) 



(27) 



(28) 



//i 



--ghost 

where we consider the change of F^ under infinitesimal coordinate transformations x 
x M + s/kV 11 . From 0(x M ) = 0(rr /M ) + 5(p and Vv( xM ) = VvC 3 '" 1 ) + ^/V> we fi n d U P to first 
order in V 1 : 



<ty. 



^ 



—a 






(29) 



and thus 



£ g host = a^^V" 01 ^D s - (D - 2)(tf 2 - a" 1 ^)^ + a" 2 **' 2 ^ V" 4 , (30) 



where V and V are the ghost and anti-ghost fields. 



The quadratic lagrangian (|2"%H251) contains mixing between the different components. 
The following field redefinition removes the mixing between ipij and ipoo and on-shell q 



5ij 

tyij = Zij + — Zqq 

"000 = Zqo 
tyOi — Zoi ■ 



£^ + Cgf + £ ghost = lx*G*"Xra + ^o;Z>Lt or ^ + VT^Y*, (32) 



The resulting quadratic lagrangian (which is still valid off shell) can be written as 

\xJ j G^ kl X kl + \zo{D^z 03 + V-Vf^V^ 
where 



Xij — Zqo 



^el r a D y5^ a D ~ 2 ^Z5 



IJ 



QijM = a Dy 6 kl p^^ a D-2^2 Z _ $ 



a D - 2 ^Z5 kl a D - 2 ^(-^Z-<!>") V 4 



D-3' 

1 .,, „,, 1 



Vll OT = a D (p. - W) (-W - -5^6 kl 

V, = a D (n s -^($)+a-VK 
Vt^ = a D ( VaP n s + 2X8l mo 



(31) 



(33) 



Note that G^ kl contains the tensor as well as the two scalar (gravitational and matter) kinetic 
operators. 



2 In the special case when D = 4, this and Eq. (140j) below agree with Ref. [2^ 
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4.1 The propagators 

We shall now construct the propagators, associated to the various modes in (1321 . In gen- 
eral this will not be possible, due to the nontrivial dependence of the propagators on the 
background fields and the mixing between the different modes. Therefore we shall restrict 
ourselves in calculating the on-shell propagators. As will be clear from the discussion in 
section 14.21 this will be sufficient to calculate the one loop effective action. It will turn out 
that kinetic operators can be written in terms of 



V„ 



-g 



D.-n d 



n 



n(n — 1) 
1 + -^ -e 



'I 



with an associated propagator 

V n iA n (x; x) = i5 D (x — x) 



(n 



)H< 



0,1,2) 



(n = 0,l,2) 



(34) 



(35) 



The operator fl34|) is however nothing but the kinetic operator for the massless scalar field, 
with conformal coupling 



£ 



n(D — n — 1 + 



n(n— 1) 



(D-l)(D-2e) 



■(1-e) 



(36) 



such that we can use the propagators calculated in the previous sections with the parameter 
is, given in ( TlBl) replaced by 



D-l 



n[D 



n 



1 + 



n(n-l) , 



2(1 - e) J 1-e 

We find for the kinetic operators for the vector and the ghost on-shell 



(37) 



Tf 3 

■^vector I on shell 



Vf° st | on shc n = ( tj^Do + ^^2?i 



(38) 



and their associated propagators: 



?'A fc 

3 vector 



-^tAj 



3 
*«^t«t = (^ iA + Og»'Al 



(39) 
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There is still mixing between zqq and 0. The on-shell part of the mixing can be removed by 
the following rotation 

X = RY 



Y 

± 7.1 



R 



A 



tan(20) 



""'.7 

X 
v 

1 

>/Acos(0) 
^sin(0) 

/ 2(D-3) 

2y/(D-3)€ 

D-3-e ' 



-V^sin^) 

^cos(0) 



arccos 



D-3 + e 



This rotation reduces the on-shell part of the term ^X GX to 



_v T r <i i kl V I 

9 ij ^diag 1 kl Ion shell 



2 lJ 



1 (±<W 



\ 



- \8^8 kl 






n 



o 









Y> 



kl 



Thus the associated propagator matrix A4diag, defined by 

^ diaa~'*idian\<L 1 <£ ) Ion shell J-0 I *£ 



T diag-/^diag\-^i •£) |on shell 



X 



is given by 



.M 



diag | on shell 



l>rs^kl 




tAr 



(40) 



(41) 



(42) 



(43) 



D-3 

To avoid confusion with the subscript diag, we have omitted the explicit Lorentz indices in 
M.diag and Gdiag- It now follows that the non-diagonal propagator matrix rs M.ui that inverts 
Qijki on-shell i s 



rs-f^Lkl on shell 



R-MdiagR | on shell 

A 2 (cos 2 (fl)A + sin 2 (#)A 2 ) 
Acos(#)sin(fl)(A - A 2 ) 



Acos(fl)sin(fl)(A -A 2 ) 
(sin 2 (£)A o + cos 2 (0)A 2 ) 



(44) 
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This finishes the construction of all the propagators and we indeed find that all modes can 
be described in terms of A n . 

4.2 One- loop effective action 

In this section we shall first sketch using a simple example how to calculate the correction 
to the Friedmann equation due to the one-loop effective action. Afterwards we shall apply 
it to the case at hand. We consider as an example a model with an action 

S = S + S x , (45) 

where x is a quantum scalar field with an action S x = S x [x] and So is the classical action 
of any background fields (including for example the Einstein-Hilbert action). For any action 
S x that is quadratic in x, 

S x = [d D x~ X T> x x, (46) 



2' 



one gets the following effective action 



i 



r=s »- il HvBipjJ =s » + 2^^)' < 47 > 

Here the trace involves tracing over the Lorentz indices and space-time integration of the 
operator at coincidence [46] and V x is the kinetic operator of the field. 

While in principle one could - at least formally - evaluate the effective action, the object 
one is eventually interested in is the effective Friedmann equation, i.e. the equations of 
motion associated with the background metric. Moreover in the present case we need to 
work under the constraint that e is constant. As long as e remains small, there is no problem 
with imposing such a constraint on the equations of motion. On the other hand, imposing this 
constraint on the level of the action typically changes the dynamics substantially. Therefore 
we shall not attempt to explicitly construct the effective action, but instead we shall directly 
calculate the effective Friedmann equation. By taking the functional derivative of the action 
with respect to the scale factor a = a(rj), we obtain the Einstein trace equation, that is the 
— (00) + (D — l)(ii) component of the Einstein equation. The second Friedmann equation 
can then always be obtained by imposing the Bianchi identity. Thus we are interested in 
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calculating 



8T 8S t 8 , . 

8a 8a 2 8a v XJ 



Va 



D-l 



D(D-1)(D-2), h2 



K 



D-l 



A + — trt?') +{D-l)p 



D 



>M 



Pm 



(48) 



+ ^ ln ^) 



where V = J d D ~ l x denotes the volume of space and we assumed that So contains the 
Einstein-Hilbert action, and matter fields with an associated total pressure and energy pm 
and pM- Notice that the quantum contribution |^-Trln (T> x ) is by definition nothing but 
Va D ~ 1 g flu (T^ l/ ), where (T^) is the one loop expectation value of the stress-energy tensor [46J. 
We now consider the calculation of this contribution. To be explicit, we shall assume that \ 
is a massless minimally coupled scalar and therefore 



V x = ^TgU 



with an associated propagator iA(x;x) which obeys 

V x iA(x; x) = i8 D (x — x) . 
Instead of considering the x field, it is convenient to use a rescaled field 



--1. 



X = a 2 X 
with an associated kinetic operator and propagator 



(49) 



(50) 



(51) 



V, 



d 2 + 



(D - 2)(D - 4) a' 2 D-2a" 



+ 



iA{x;x) = a^ a? iA(x;x) 



(52) 



We shall see explicitly later that - up to a D dimensional divergent delta function that does 
not contribute in dimensional regularization - we have that 



ATYln(Z> x ) = ATYl, l( P x ). 



(53) 



We shall now show how to calculate the trace logarithm contribution to ff48l due to the field 
X- To be precise we give the exact coordinate dependence of each term indicated by x^, y* 1 , 
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and z^ 



2 5a(z< 



-Trln(P x (x)) = lp(h{ x - y )(*V x {x) 



- J d D x J d D yA(x;y)S D (x-y) 



(D-2)(D-4)a'{x ) 2 D-2a'\x 



,"(„o-^ 



+ 



l -v 

2 



2 a(x ) 3 

(D-2)(D-A) a'(x ^ 



a(x 



0A2 



5(x° - z°) 



D 



2 a. 

{D-2)(D-A)a'(z ) 2 
2 a{z°f 



M° 9 *" Hx0 ~ z ° H ^) a '' Ax °- z0) 



2a(x° 
D-2a"(z° 



2 a{z u ) z 



(D-2)(D-4) a /a'(z°) 



2 
D-2 



d 7 o 



m iAiz :) 



+ ^A 



z v a(z° 



-iA(z u :z 



0. „(h 



Va D - l U z iMz\z) 



(54) 



where in going from the second to the third step we used the delta functions to change d x o 
to —d z o and we used that the propagators at coincidence are a function of time only. In the 
last step we used (I52p to rewrite iA in terms of iA. 

Since the trace of the one loop expectation value of the stress energy tensor for a scalar 
field is given by 



(0\T%\0) 

we see that indeed we have that 

i 5 



2- D 



(D- l)nniA(x;x), 



--Trln(V x ) = Va D - l g^(T llu ). 



(55) 



(56) 



Comparing this with Eq. (Jl"8"j) we see that this is exactly what one expects, justifying thus 
our procedure. The above example shows that indeed one can use the functional derivatives 
with respect to a to evaluate the quantity that interests us. It also shows that the rescaling 
(I5T]) does not influence the final result in the context of dimensional regularization. Notice 
that (1481) is just an equation of motion. Therefore, after the variation is performed, in 
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the second line of (1541) . we can safely evaluate all the quantities appearing on-shell. Hence 
we only need the propagators iA(y; z) on-shell. This justifies our on-shell diagonalization 
procedure, based on which we constructed the propagators. 

Now we apply this technique to the case at hand. The only difference is that there are more 
quantum fields. This however does not change the procedure. Up to an irrelevant constant, 
the effective action can be obtained from 



exp[ir] = / (Vh^) (V<j>) (VU a ) (VU a ) exp (i(S {0) + S (2) ) 



(57) 



where U and U denote the unrescaled ghost fields associated with the graviton field h^ u . 
When written in terms of the rescaled fields this can be recast as lj 



exp[*r] = / (Vz tJ ) (Vz 0i ) {Vz 00 ) {V<j>) (VV a ) (VV a ) exp (i(S^ + S&) 



exp 1 1 



: S (0) 



Tjghost 
U af3 



(5? 



detpl^) det(G^) 



where 



S-(°) = J d D x^(-(R -(D- 2)A) - i(9$) 2 - V($) 
5( 2 )= fd D x£( 2 \ 
and £ (2) is given in (JB2J) and (<9$) 2 = -$ /2 . From Eq. (jHHJ) we obtain 



(59) 



r = S<® + l -Tr ln[Pl tor ] + l -Tr HG^ 1 ) - *Tr \n[V a0 
= S {0) + T 1L , 
from which we find the expression equivalent to (1481) to be 

sr 



hosti 



(60) 



_ 5So STil 

5a 5a 5a 



Va 



D-l 



D(D-l)(D-2) , ff2 



K 



_L_ A + ^a- l H') + (D-l)p M - p M 



5T 



XL 



5a 



(61) 



3 Our field redefinition ([31]) has a Jacobian equal to one. Furthermore, the rescaling by a 2 of rp with 
respect to the 'true' graviton will contribute to the effective action as a Z?-dimensional delta function S D (Q). 
Such a term does not contribute in dimensional regularization. 
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where pm and pu are the pressure and energy density associated to the background scalar 
field matter, given by p M = \& + V($) and p M = \& - V($). 
The one loop contribution can be written analogously to (j54p as 

1L = d D x ( -\A vcctoI ](x-x)—V ij + -[■ Mi,,]— G ijkl -i\ a &? ](x-x)—V ghost \ 

(62) 
The functional derivatives should naturally be taken on the off-shell kinetic operators. How- 
ever, as soon as these derivatives are taken, we are simply left with an equation of motion. 
It is therefore completely valid to impose the background equations of motion at that point. 
Therefore the propagators in (162]) can be evaluated on-shell and thus we can use the propa- 
gators as they are calculated in the previous section. 
Instead of using the kinetic operators as given in (13"3"j) . we rescale all of our fields as 

V-a 1_iV2 <V ; 0-a 1 "^, (63) 

which is identical to the rescaling in Eq. (15T1) . This rescaling changes the effective action 
(I58p by a D dimensional coincident delta function that does not contribute in dimensional 
regularization. With the following identity 



-<?a 



V a "d a d, + ^^-e)a^ 



(64) 



we can easily calculate the kinetic operators of the rescaled fields. We will indicate these 
rescaled kinetic operators with a hat. The associated propagators are also easily obtained: 

A(x; x) = (aa) x ~ D/2 A(x] x) . (65) 

We now shall consider one functional derivative in detail. The others are calculated similarly. 
We follow the same procedure as in the example considered at the start of this section. Our 
calculation proceeds analogously as in Eq. (1541) . We once again insert the explicit arguments 
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i] = x° and fj for the two coordinates respectively 

a X ijl\kl\X]X) , . -^tensor W J 

= Jd D x lJ A kl (x;x) J ^ ) L a "d a d, + ^(y " e(v))a(v) 2 H(v) 2 - a{f,fW 
V2 4 / 



= / A ( ,A H (.;^- __( 2(3D _16)( -|1- - JL) 

+ 6(- ^ a( '' ) ") S ) 
0(77) ((9r]) 2 a(fj) 2 J 

+ 2(D - 2)Aa(r/) + 2/d/($)a(?7) ) 0-5 lk 5 jl - U ij 8 k{ \8{f, - rj) 



2 0(77) (<9r/) 2 
+ ±(£> - 2)(3L> - 10) (H A + (1 _ e)H(f,ya(v\ 

+ 2p - 2)Aa(f,) + 2KV($)a(fj) j (^ fc <$* - ^<5 fc ^A H (x; x) 



3(0-2)1 1 

a " 2 



V( -^-^-^ 2 H--(L>-2)(3L>-10)^ 



(66) 



+ -(D - 2) 2 (7 - 3e)tf 2 a j (-,W - -^5 kl ) tJ A kl (x;x) . 
In step three we integrated by parts and dropped the boundary terms because of the delta 



function. In step four we used the background equations of motion (j2"0|) for $. This is justi- 
fied, since corrections to those equations will be of order one loop and the above expression 
is already at order one loop. Therefore the error one is making is of two loop order. 
The other functional derivatives we need are calculated similarly 

' f A^Ux; x)^Vl ctOT = ( ^\ 2 -\{D- 2) (3D - 2)Hd v 



T r I ' — JjVuuloi \.~; -~/ r — vector I r, 

V J oa \ 2 a 

* j, vector \%i % ) 



- i(£> - 2) ((3D +2)(1 - e) + 4(D - 2))aff 2 j^A, 



1 'jCW; ^^Acaia, = ^r4r ^k + l( D " ~ 4 )^ 



+ i(D - 2) ((7D + 2) - (3D + 10)e)if 2 a )>f (1)1) (x; x) 

.M (2 , 2) (x; x)^, = f ^^ - ^(D - 2fHd, - (±(D - 2) 2 (1 - e) 
+ 4(D-l-e)e)aiJ 2 J>l (2i2) (x;x) 
1 ^Af host *-V§* = (^\dl - \{D - 2fHd v - \{D - 2) 2 (1 - e)H*a) Va ^ 



(67) 



7 7 ^ghost 5a -a/3 "V 2 a" 2 V ; " 2 V > y ' J lap ghost 

+ ((D - 2)-^ 2 + 4(D - 2)Hd v + 4(D - 2)(1 - e)if 2 a)<fc a Af host 

= ^f_~ 2)e ((D - 2)#d„ + (3D - 4 - De)H 2 a^M (1 , 2) (x; x) . 

Here we indicated with A^( n , m )(a;;x) the (n,m) component of the propagator matrix .M 
( |44|) . The last thing we need before we can calculate the one loop contribution (|62|) are the 
propagators at coincidence and their derivatives. Since all propagators are related to the 
propagator A n , we only need that one. 

4.2.1 Infinite volume contribution 

The propagators that we use are split in two parts ffl4|) : the infinite volume part and the 
corrections due to the infrared cut-off. 
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We first consider the infinite volume contribution HT5]) . Taking the y — > limit of (j!5p and 
dropping the D dependent powers of y that do not contribute in dimensional regularization 
we obtain (we drop the subscript oo) 



iA n (x;x)=a D 2 iA n (x; x) = |1 — e 



D-2/ tt\D-2 



aH) 1 



r(i-%)r(^ + un n )r(^ 



-vd. 



(4tt)t 



T(\ + v Dyn )T(\ -v D>1 



iA n (x; x) = Ha(D - 2)(1 - e)zA„(x; x) 



—J *A n (x; x) = H 2 a\D -!)(£>- 2)(1 - e) 2 *A n (x; x). 



(68) 



We can now collect all the terms of (J52J. Using (pj) . (l44|) . IjBSj) . (IB7JI and (l68l) we obtain for 
the vector and the ghost 



_ L.A vector (r- r)—T> ij 
V / 2 l j {X,X) 5a vector 



(D - \){D - 2) 2(D - 1)(D + 2) 
(D 2 + D + 2)e - (D - 1)(D - 2)e 2 



aH 2 i Ai(x; x) 



1 
V 



-^ Af host (x; x)^P$° st =^ 2 a^ - 1) 2 (D - 2) 2 e(l - e)zA (x; x) 



- -(D - 1)(D - 2)(1 - e) h(D + 2) - 3(D - 2)e)aiJ 2 zA 1 ( 



x; x) 



(69) 



The last term evaluates to 

\ J ^jM kl (x; x)y a G ijU = i I iyAuCar; aO^EgL 

+ -A 2 ( COs(#) 2 «A + sin(#) 2 *A 2 ) ^scalar 

+ Acos(#) sin(^)(^A - %A 2 )— ( yfc- — -Z - ^") 

oa\ D — 6 

+ i(sin(^A + cos(#) 2 *A 2 )^. 
The contribution from the tensor is 



(70) 



-D(D - 1){D - 2) 2 (1 + 3D) - 3(D - l)e 



eaH iAn(x:x) 



(71) 
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while the terms in ( 1701) multiplying iA contribute as 



A(D - 3 + e) 



(D -!)(£>- 3)(£> 2 - 8D + 4) + [D(D - 2)(D 2 - 11 D + 14) + 8]e 



+ (D - 1)(D - 2){D + 2)e 2 



a/7 zAnfx; x) 



(72) 



Finally the terms in (ITUj) multiplying «A 2 contribute as 
1 
4(L> - 3 + e) 

D 4 - 5D 3 - AD 2 + 24L> - 32)e 2 4 ( D 6 - 5D' 2 + 8£> + 4 )e a 



4(Z> - 1)(D - 2)(Z> - 3)(D + 3) - (W 4 - 20L> 3 - 9L> 2 + 68£> - 36)e 



(73) 



aH 2 iA 2 (x; x) 



Putting (J69H73J) and ( 1621) together and expanding the result around D = 4 we obtain the 
following non renormalized one loop contributions to the Friedmann trace equation, 



1 5T° 1L e(Q3e 2 + 2e-W5) n , 2 „ 2 ,J 2 . /(l-e) 2 # 2 

K^rf = 64^(1 + e) (1 " e) <** " 1} (f! + 1E + ln ( 4^~ 

+ 327r2 " e)2 (l " e) 2 (4^ 2 - 1) (93e 3 + 906 2 - 169e - 122 



1 5T\ L 3(9e 2 - 7e - 6) 



Va D ~ x H^ 5a 



647T 5 



(l_ c )3(4^_i)(_z_ + - L + b. 

4z/i^ - 1 



2tt2 



47T/i 2 



+ ^(tt + ^i) + ^(tt - ^i) + 4 a , 



Auf - 1 



+ ^(l- e ) 2 (4^ 2 -l)(51 e 2 -17 e -54 



2 / (1 — f) 2 H 2 



i ^ = (5e-6)(e 2 -2e-7) 



4 X /I X 4Z/ 2 ^n - 1 

+ <Mo + "*) + ^(o - "2) + 4 ^7T 



4^ - 1 



+ — -[— -(1 - e) 2 (4z, 2 - 1) ( 3 e 4 + 13e 3 - 83e 2 + 35e + 3o) . 



(74) 
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Here is n indicates the index is n given in (157)1 with D = 4 and v' n is -^vr>,n 
Euler constant and we used the following expansions of the propagators (168 



D=4 



. 7£ is the 



a/7 iA„(x: x) 



rei - -' 



o^ll-e 



rfi - ^ 

D-2ttD l \ l :> 



(47T 



2 



L>-3\ 2 



x 



1 + £ y i (^(^ + z/ ") + H^-^ 



(75) 



2rr4 



e W 



(4tt)- 



x 



167T 2 



In 



47T 



(76) 



plus terms that vanish in D = 4. Here ^(z) = (d/dz) ln(r(^)) is the digamma function and [i 
is an arbitrary renormalization scale introduced for later convenience. If we use the explicit 
expression for v n , we can add all terms together and obtain for the functional derivative of 
the effective action 



ST 



\L 



^\l 



+ 



W\l 



+ 



5T 2 1L 



Va D ~ l H A 5a Va D - x H A 5a VaP~ x R^ 5a Va D ~ l H A 5a 
e(l86 - 149e - lie 2 + 10e 3 ) ^d-a 



8vr 2 



D-A 



167T 2 



108 + 62e - 153e 2 + 27e 3 



186-149e-lle 2 + 10e 3 



x l-f E + ln(- 



2tt2 



1-efH 

Att/j 2 



e(63e 2 + 26-105),, s2fA 2 „J , ,1 , ,,1 



"a) 



(77) 



3(9e 2 -7e-6) ., .„., 2 1X / ,,1 ,.1 
— Hi " e) 2 (4z/ x 2 - 1) ( ^(- + I/O + ^(- - //., ) 



64tt 2 

(5e-6)(e 2 -2e-7) 
64tt 2 (1 + e) 



(l- e ) 2 (4z, 2 -l)(^ + z, 2 )+^-z/ 2 ) 



We kept the digamma functions in terms of is, since then it will be more clear how to add 
the correction terms to regulate the infrared. In ( 1771) . all infrared power-law divergences 
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are (incorrectly!) subtracted by the automatic subtraction of dimensional regularization. 
The infrared logarithmic divergences are still there and they appear through the poles in 
the digamma functions. These issues will be corrected by adding the correction terms, but 
before adding the correction terms, we shall first renormalize the expression (177j) . 

4.2.2 Renormalization 

The contribution (177j) contains a l/(D—4) divergence and therefore needs to be renormalized. 

If we take our approximation that e = constant literally, the divergence is of the form constant 

xif 4 a Z3 ~ 1 . However in a more realistic treatment, e is a dynamical parameter and our result 

is expected to be correct up to zeroth order in e. Such an approach, using a space-time 

of locally constant e is a generalization of the often considered locally de Sitter space-time 

[39J [50]. In this more realistic case, the e structure of the divergent term should be taken into 

account and subtracted accordingly. Therefore if e is varying slowly enough, such that (1771) 

remains approximately valid, we still need a counter lagrangian that produces the same e 

structure as the divergence in (I77j) . in order for the theory to be renormalized at all times. For 

this purpose many terms can be used, and indeed many terms are reported in the literature 

[TJ1 [9], H7J . However we cannot simply use these terms, since counterterms are dependent on 

the gauge fixing used [48J and as far as we know, no general calculations have been done 

using our gauge fixing (j2"7|) . Calculations using (1271) have been done however in the special 

case of de Sitter space [2H EH) considering both scalar and graviton loops. From these 

works and also e.g. from [13] it follows that the one loop contribution in this limit should 

be finite. Since the de Sitter limit is e — > 0, this agrees with (!77|) . To ensure that the one 

loop contribution due to gravitons in the de Sitter case vanishes we need in that limit a 

counterterm 

A 



-»[*'- ij-i.) (* 



which in our more general case becomes 

Ccti = V^a (R- -^ (kV{$) + (D- 2)a)) ' , (79) 



D 



where a® is a constant. This follows from the fact that the cosmological constant can always 
be seen as a part of the scalar potential and thus for an invariant counterterm they should 
come together. Moreover from [H] it follows that we also need a counterterm of the form 
V —gH 2 ng^ v '(dfj,(j))(d u (j)) to vanish. From Eq. ([9]) it follows that for our case this generalizes 
to 

C C T2 = V^ga^Rg^ - DR^) K {d^){d v $) , (80) 



23 



with ax a constant. Finally, also from [44] it follows that a counterterm y/— g/t(n$) 2 should 
not appear. 

Therefore a reasonable choice for the counter-lagrangian is 



£ c = \F I g 



m>(r - q^( kV + ( d ~ 2 ) a )) 2 + «i (^ - ^i?^)«a^^$ 



<9$ 2 ' ° * ^'^"-' ^ 



We stress that the counter-lagrangian (1811) /or t/ie purpose of this calculation could be chosen 
differently. There are many other terms with the correct dimensionality that could have been 
used [1?]. Since the divergence (ITT)) gives only four constraints (one for each power of e), 
we at present can fix only 4 coefficients. This does not mean that the counter-lagrangian 
is arbitrary. Different types of calculations could fix the counter-lagrangian uniquely, as 
it is for example done in Ref. [9]. However apart from the two cases mentioned above, 
we do not know of any calculation in our gauge, which we could use to further specify 
our counter-lagrangian. Thus the 'true' counter-lagrangian corresponding to the theory will 
probably contain different counterterms than (1HT1) . However, these different counterterms 
do not change the conclusions of this paper. The only effect would be that in Eq. (|85l) the 
origin of the /3j's changes, but not the fact that they are essentially arbitrary. The terms in 
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our counter-lagrangian (1ST]) contribute as follows to the Friedmann trace equation 
d D x^g(R - 



Via 



D 



D N 2 

' -(kV+(D-2)A) 



,D-1 rr4 



D-lrj4 



# 4 (D - 2)e(8e(2 + 3e) + D 2 (A + 9e) - 2D(2 + 13e + 12e 2 )) + C(e') 
Wl6e(6 + 7e - 9e 2 ) + 4e(26 + 37e - 30e 2 )(L> - 4) + 0((D - 4) 2 , e')) 



1A 

Via 



d u xK^j(Rg^ - DR» v )(d^)(d u §) 



D-IttA, 



1 8 

Via 



lD 



d u x^g 



H\D - 1){D - 2) 2 e((D - 1)(D - 6e) + 6e 2 ) + £>(e') 

rr ^ Wl44e(l - e)(2 - e) + 24e(23 - 30e + 8e 2 )(D - 4) + (9((D - 4) 2 , e 
9 2 ^($) 



2a 

„D-lrr4 



a 



^$2 
D-l rr4 



-R 



n-izjA 



H 4 (D~l)(D-l-e)2e(D(D-2) - 2(3D-4)e + 12e 5 
W24e(3-e)(2-4e + 3e 2 ) 
+ 4e(51 - 88e + 53e 2 - 6e 3 ) (D - 4)) + o((D-A) 2 , A 



0(t 



]_6_ 
VTa 



-— / d n xKV=gg^(d^)(d 1/ $) 



d 2 V($) 



Aa D ~ 1 H\D-2) 2 e 2 (D-l-e) + O(e') 



(82) 



13) 



,D-1 rr4 



W - 16e 2 (3-e) - 16e 2 (4-e)(D-4)) + o({D-A) 



V , 



where we used ^ and once again we used the background equations of motion f[2"Uj) and (12TJ) . 
We find that all divergencies cancel if 



a 



37 /j, 



D-4 



9Q0tt 2 D-A 

D-4 



a 2 



fJ 



288tt 2 D-4 



+ a f ; 



+ 4 



a i 



a 3 



49 fi 



D-4 



64Qir 2 D-4 



+ a{ 



43 /i 



D-4 



54) 



480vr 2 D-4 



+ 4, 



where the af (i = 0, 1, 2, 3) indicates a possible finite part. Adding the contribution from the 
counterterms ( l82l) to the one loop contribution ( 1771) . we obtain the following renormalized 
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contribution 
a 3 V 5a 



H 4 

167T 2 



P 1 e + p 2 e 2 + p 3 e 3 + p 4 e 4 



e( 186 - 149e - lie 2 + 10e 3 ) (' In 



1-eYH 



2tt2 



A'KjJ? 

e(63e 2 + 2e-105). i . 2 , A 2 1N /.,1 , ,1 

+ 4(1 + e ) < l - e ) ( 4z7 o " 1) ( ^ + Vo) + ^ ( 2 " Uo) 

+ - ' 2 ~ 4 Te " 6) (1 " e)2(4 ^ " 1} 64 + Ul) + ^\ ~ Ul) 



15) 



3(9e 2 - 76 - 6) 

4 l 

(5e-6)(e 2 -2e-7) 



(1 - e) 2 (4^ 2 - 1) ^ + ^) + ^(o - ^) 



4(1 + e) J \ '2 " '2 

Here the parameters /3, (i = 1, 2, 3, 4) are given in terms of the finite coefficients oq 



Pi = 16vr 2 x 12 



8a f Q + 24a{ + 12a^ 



186 7£ ; + 



1727 



02 



16tt 2 x 4 



28a f - 108a{ - 84a{ - 12a{ 



p 3 = 16vr 2 x 4 



- 36a f + 36a{ + 78a f 2 + 4a{ 



+ 1497s - 
+ 117s + 



5969 

~9~ 

10457 

45 



(86) 



/3 4 = 16vr 2 



-72a: 



-10 Ti 
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All /?j's (i = 1,2,3,4) in Eq. (1851) are free parameters that remain undetermined until they 
are fixed by experiment. 

4.2.3 Correction terms 

The effective action (1851) is divergent for half integer values of the parameters z/j > 3/2 (notice 
that the pole at z/j = 1/2 is cancelled due to the pre-factor) 



Ui = 3/2 + N 



N = {0,1,2...}. 



(87) 



The reason is that the propagators used do not describe the infrared physics correctly and 
we need to add the correction terms 5iA^ and 5iA N as given in ffT6|) . For concreteness we 
shall consider the late time behavior in an accelerating space-time, since in that case we 
expect the most significant backreaction. In an accelerating space-time, 77 goes to zero at 
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late times, and thus we do not care about the SiA N corrections, since they quickly become 
insignificant in that case. Since the correction terms are ultraviolet finite, we can put D = 4 
in all terms. Following the same procedure as in section I4.2.1[ we find for the coincident 
limit and its derivatives 

8iA N , n (x;x) = A N , n (aH) 2 z 2 N+3 - 2u " 

d ~ * 

— 5iA Ntn (x;x) = -(1 + 2N - 2v n )(l - e)(aH)SiA N>n (x; x) / gg N 

—J SiA N>n {x; x) = (l + 2N- 2u n ){2N - 2i/„)(l - e){aH)8iA Ntn {x\ x) , 

where v n is given in ( l3"Tj) and a subscript n implies that the quantity is evaluated with v = v n . 
We also defined 

A - 1 l T(v n -N)T(2v n -N) 2 

A "' n " 4vr5/2 3 + 2iV - 2^ n T(i + v n - N)T(N + 1) l j ' l j 



:| -4(3- e ) e + (l-e) 2 (l + 3 e )(2^ + 2z/ 2 )-(3 + e )(l- e )(l-3e) 

+ JV(1 - e) (3 - 4z/ - 8e(l + u ) - 3e 2 (l - 4i/ )) ] aH 2 5iA Nfi (x; x). 
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(90) 



We find for the corrections due to the vector and the ghost (jo 
^ JfaAf^Mx; x^vLtor = 34 + 7(1- c) 

- 5(1 - e)(l + 2N - 2v x ) - (1 - e) 2 {N - v t ){l + 2N - 2v x ) J aH 2 5iA NA {x; x) 

\\~ ^ a K^ s Mx; x)yj>tt = 6(1 - e)(JV - 1*) 
x ((2z/ - 2iV)(l - e) + e - 3)aH 2 5iA N , (x; x) 
+ 6(1 -e)(N- ui) ({2v! - 2N)(1 - e) + e + l^ofl^iA^ar; x). 
The tensor contribution to flTDl is 

- 6(3(i/o - N)(l - e) + 2) (2(i/„ - iV)(l - e) - 3)aH 2 «JzAjv >0 (a;; a) , (91) 

while the terms in (1701) multiplying 5iAjv,o contribute as 
1 



(92) 



Finally the terms in (17DJ1 multiplying 5iA N:2 contribute as 

-^ I (1 - e) 2 (3 + e)(2N 2 + 2v\) + 2(3 + 5e(l - e) + e 3 ) 

+ (1 - e)(3 + e) ((5 + e){l - N) - 4(1 - e)Nv^ J aH 2 SiA N>2 (x; x). 
We add the corrections together to obtain the following three contributions 



1 6lf_( e (63^ + 2e--105) + l-e 



(93) 



VaP^E^ 6a \ 1 + e 1 + e 

5iA Nfi 
H 2 



x ((u -N)(l- e)(23 + 21e) + (4 + 3e)(3 - 7e) 



<5r (1) , 

v ' 3(9e 2 + 7e + 6) 



VaP~ x H^ 5a 

+ 9(l-6)(3 + 2iV-2r, 1 )((^-iV)(l-6)-e)j^fi 

1 oTg } /(5e-6)(e 2 -2e-7) 



W^tf 4 <5a \ 1 + e 

- <I^±i>(3 + 2N - 2„ 2 )((„ 2 - JV)(1 - e) + 4 - e j\ **£*, 

(94) 

where we have written £TJy for the contribution to 5T N that multiplies Ajv iT1 . 

From these three contributions we see explicitly that indeed when v = N + 3/2, the 
corrections have the correct prefactor to add up correctly to cancel the divergence in the 
digamma functions (j85|) . Inserting (|85|) and (1941) into floTT) at D = 4 we obtain the final one 
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loop corrected Friedmann trace equation 

^ ((1 - l -e)H 2 - 1a) + 3 PM - PM + -^ J /3 l£ + (3 2 e 2 + /3 3 e 3 + &e 4 

( 2 ,\, /(l-e) 2 # 2 \ 
- ef 186 - 149e - lie 2 + 10e 3 J In ( K \ J 



e(63e 2 + 2e - 105) , x9 

+ — -(1 - e) 

4(1 + e) l J 



x ((H 2 - 1) (^ + H>) + ^ " -o)) - £ (g^/ i* 



, 3(9e 2 -76-6) 



x ((4*? - 1) (^ + v.) + ^ - ^i)) - £ (ff^/ ^w) 



( 5e _ 6 )( e 2_ 2e _7) 



4(1 + e) 
x ((^ 2 2 - 1) (i>{\ + m) + ^ - V2)) ~ £ ( ff ( 6 i 4 - e ))2 ^^ 

+ £(3 + 2iV - 2i/„) f ^ ((i/o - iV)(l - e)(23 + 21e) + (4 + 3e)(3 - 7e)) J # 2 <5iA,v,o 
+ £(3 + 2iV - 2 ^i) ( 9(1 - e ) ((*"i - #)(1 - e ) " e ) ) H 2 5tA NA 

- £(3 + 2iV - 2z/ 2 ) f (1 "^^ e) ((^2 ~ N)(l - e) + 4 - e) J H 2 6iA N , 2 = . 

(95) 

This equation is our final result of this section. It presents the one loop quantum corrected 
Friedmann trace equation, in the presence of both graviton and scalar fluctuations. 

4.3 Discussion 

Having found the one loop corrected Friedmann trace equation, we can ask the question 
whether the quantum effects calculated will have any significant effect. We immediately 
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see from (1951) that the quantum contribution is suppressed by a factor H 2 /m 2 . Because of 
this suppression, the one loop contribution can only become relevant for the dynamics of 
the Universe if there is a significant enhancement. In other words: for any significant effect 
the one loop correction terms should grow in time with respect to the classical tree level 
terms. In an accelerating universe zq goes to zero at late times. From flHHj) we see that the 
corrections <5iAjv,n are proportional to H A z^ +3 ~ 2Un and thus this contribution might grow in 
time with respect to the background (which scales as H 2 ), making the initially tiny quantum 
effects significant at late times. The question whether this growth is significant is equivalent 
to the question whether the quantity 

H 2 z 3 + 2N-2» (96) 

grows in time. In an accelerating universe the fastest growing term has N = 0. Since 
zq = —krj this term grows with rj as 

oc ip^~ 2v . (97) 

Thus we see that for the three different z/ n 's contributing to ( |95l) this becomes 

oc r]^~ 2uo = r]° 

oc r^- 2 " 1 = V 2 (98) 



3-e 
ll-e" 



2^2 _ ^4 



OC I] 1 -* ' = 7] 

These general expressions fail near the N = pole of the digamma functions (meaning 
Pi = 3/2). Here one will typically pick up an additional logarithm, ln(z ) from the correction 
terms. Thus we see that the uq contribution, in the presence of such an additional logarithm, 
actually grows in time. However u = 3/2 implies e = and we see from (I9"oT) that this 
contribution is cancelled by the pre-factor. Thus none of the contributions actually grows in 
time. 

The fact that the pre-factor is such that near de Sitter space the logarithmic growth drops 
out, is also in concordance with results obtained in [23J. This property however does not 
have to stay true in general. For example if higher order loop corrections are taken into 
account, this structure might change such that while the correction terms are extremely 
small, their contribution grows in time, making them significant at late enough times. Also 
when considering a different background geometry, one might get growing effects even at one 
loop order. In [161 02] f° r example one loop contributions have been considered in chaotic 
inflationary models. Although those works found a secular growth, in [21] it was concluded 
that this growth disappears if one considers truly gauge invariant quantities. 
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5 Conclusion 

In this paper we have calculated the one loop contribution to the effective Friedmann equa- 
tions due to graviton and matter quantum fluctuations. The background space-time we 
use is FLRW with the additional constraint that e is constant. In such a calculation many 
delicate issues have to be taken into account. First of all there is mixing between the gravi- 
tational and matter degrees of freedom. This mixing can (at least on-shell) be removed by 
field redefinitions (131]) and fj40|) . Secondly, there is the issue of gauge freedom, which we can 
deal with in the standard way by adding a gauge fixing term (1271) and the associated Fadeev- 
Popov ghost lagrangian (1301) . Finally, to calculate the one loop contribution, one needs to 
know the propagators. The propagators needed will however typically be infrared divergent. 
To regulate this infrared divergence we assumed that the spatial sections of the universe are 
compact. This effectively implies that the propagators should be calculated using an infrared 
cut-off k . Using this approach the scalar propagator was constructed in [I]. In the present 
work we have constructed the graviton and matter propagators by making use of the same 
method. All these technical issues can thus be resolved, leading to our final answer for the 
one loop Friedmann trace equation (j95p . Although the final answer might look intimidating, 
we can essentially write it as two types of contributions. One is H A times a constant and the 
other is H A times some power of z$ = \rj\ko- Since the tree level contribution is of the order 
of H 2 i7ip, we see that - as expected - quantum corrections to the Friedmann equation are 
suppressed as H 2 /ml. There are thus only two possibilities for the quantum corrections to 
have a significant effect. Either the coefficient in front of H 4 is large, or any of the powers 
of zq grows fast enough to compensate the suppression after a sufficient amount of time. 
Looking at equation (|95|) we see that for typical values of the parameters, the coefficient will 
be of order one. One might think that the digamma functions might grow large. Indeed, 
these functions have a simple pole when the argument is a negative integer. However, as is 
shown in [TJ, these poles are due to a logarithmic infrared divergence. These poles should 
therefore be subtracted by the correction terms, and in [1] it is shown that this indeed is 
what happens. Thus the only possibility for a significant backreaction of quantum fluctua- 
tions on the background space-time is if the correction terms grow fast enough in time. In 
section 14.31 we showed however that this will not happen. The only possibility where the 
powers of z grow fast enough is the special case of de Sitter (e = 0). However in this case, 
the growing power of z drops out due to the vanishing of its prefactor. 

It is however not clear whether these conclusions will also hold in more general cases. 
Generalizations could include for example studying a more general class of backgrounds, 
using different regularization schemes for the infrared. For example in [45J a scheme using 
mode matching is introduced. This scheme gives similar results for accelerating space-times 
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as the scheme used here, but the results differ for decelerating space-times. Another general- 
ization would be to include two or higher loop contributions into account. Indeed studies in 
de Sitter space show, for example, that at two loop order the backreaction due to gravitons 
might indeed grow in time [SI HSl IE]- This paper should be seen as a first step in trying to 
understand backreaction on a more general background space-time. 
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